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The natural notion of almost perfect state transfer (APST) is examined. It is applied to the 
modelling of efficient quantum wires with the help of XX spin chains. It is shown that APST occurs 
in mirror-symmetric systems, when the 1-excitation energies of the chains are linearly independent 
over rational numbers. This result is obtained as a corollary of the Kronecker theorem in Diophantine 
approximation. APST happens under much less restrictive conditions than perfect state transfer 
(PST) and moreover accommodates the unavoidable imperfections. Some examples are discussed. 



1. INTRODUCTION 

The design of models to transfer quantum states be- 
tween distant locations is of relevance to many quantum 
information protocols. In recent years, spin chains have 
been proposed [3], [1] as possible basic systems in the 
construction of such quantum channels. In these devices, 
the state of the qubit at one end of the chain is transfered 
to the qubit at the other end after some time. A key ad- 
vantage of this method is that it minimizes the need for 
external interventions as the transfer is realized through 
the intrinsic dynamics of the chain. 

The efficiency and reliability of quantum wires must be 
high. Ideally, it is wished that the probability of finding 
the initial state as the output be one; when this is so, one 
speaks of perfect state transfer (PST). It has been shown, 
in particular, that PST can be achieved in spin chains 
by properly engineering and modulating the couplings 
between the sites [1], [11]. 

In this context, a question of practical importance is 
that of the robustness of these ideal transfer properties in 
view of the unavoidable manufacturing and experimental 
deviations from the theoretical specifications. 

There are many sources of errors: nonsynchronous 
or imperfect input and read-out operations, fabrication 
defects, additional interactions, systematic biases etc. 
What their influences on the fidelity of state transfer are 
[16], [22] and how to correct or circumvent them [6] , 
[20] has been the object of various studies. We shall here 
relate particularly to the errors imputable to the quan- 
tum wires and measurements. In [9] (see also [7], [8]) the 
imperfections in the production of the device are mod- 
eled by adding random perturbations to the couplings 
and magnetic fields of a chain with PST. In [4] (which 
uses methods similar to [17] and [23]), modulated chains 
are coupled to boundary states to make transfer more ro- 
bust against imperfections which are also randomly sim- 
ulated. We shall adopt in the following a complementary 
approach. 

We take for given that the manufacturing of the chains 
will not be perfect and that the precise PST require- 
ments will not be met. Notionally, we assume that the 



imperfections are static. We then consider under which 
conditions will state transfer be almost perfect. In other 
words, instead of examining the effect of perturbations 
on PST chains, we shall readily attempt to characterize 
the chains that somehow incorporate defects and come 
"very" close to achieving PST. 

The simplest spin chains exhibiting PST are governed 
by XX Hamiltonians with nearest-neighbor interactions. 
We shall confine ourselves to these systems in the fol- 
lowing. Given that the total spin projection commute 
with the Hamiltonians, many of the state transfer prop- 
erties can be obtained by focusing on the 1-excitation 
sector of the state space. The corresponding restrictions 
of the Hamiltonians are tri-diagonal Jacobi matrices J 
whose entries are the couplings and magnetic fields of the 
chains. Naturally, these J are diagonalized by orthogo- 
nal polynomials (OPs). The perfect transfer of a single 
spin up from one end of the chain to the other is possible 
only if J is mirror symmetric; moreover it puts strong 
requirements on its spectrum. The theory of orthogonal 
polynomials is very useful to obtain these results. 

The same framework will be adopted to determine the 
conditions for almost-perfect state transfer APST. The 
question now is: under what circumstances are there 
times for which the transition probabilities from one site 
to another can be made to approach 1 . Mirror symmetry 
will again be necessary as we shall see. It will be suf- 
ficient however for the spectrum of a mirror-symmetric 
J to obey conditions much milder than in the PST case. 
Interestingly, these results will follow from the Kronecker 
theorem in Diophantine approximation. 

The remainder of the paper will proceed as follows. 
We shall first review in Section 2 how 1-excitation state 
transfer is realized in XX spin chains and how the eigen- 
states are obtained with the help of OPs. Almost perfect 
state transfer will be defined in Section 3 and necessary 
and sufficient conditions for its occurrence will be deter- 
mined. In Section 4, it will be shown how various chains 
with APST can be obtained from a parent chain with 
APST through a procedure referred to as spectral surgery 
[23]. Section 5 offers a number of illustrative examples: 
situations where APST occurs and cases where neither 
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PST nor APST are possible. In a significant model, the 
times for which APST is attained can be estimated show- 
ing in this instance that a good approximation to PST 
is obtained in finite time independently of the size of the 
chain. 

Let us finally draw attention to recent related and com- 
plementary publications that have made use of the same 
tools from number theory in the case of spin chains with 
uniform couplings. In his thesis [5], using the Kroncckcr 
theorem, Burgarth has shown that the Heisenberg XX 
spin chains with prime lengths have arbitrary good fi- 
delity. In [12], Godsil introduced the notion of pretty 
good state transfer on graphs which is equivalent to our 
definition of APST. Our terminology is dictated by the 
fact that in these instances the transition amplitudes are 
almost periodic functions. Recently, the characterization 
of the uniform XX Heisenberg chains which admit this 
pretty good state transfer was analyzed using again num- 
ber theoretic methods [13]. 

In dealing with the Heisenberg XX spin chain with 
uniform coupling and zero magnetic fields, the mirror- 
symmetry of the one-excitation Hamiltonian is de facto 
ensured as this operator coincides with the adjacency ma- 
trix of the (N + l)-path. As will be seen, this symme- 
try requirement comes into play in relation with APST, 
when more general couplings and magnetic fields are con- 
sidered. It is remarked in [12] that in order to get a good 
approximation to perfect state transfer, the waiting time 
must be very large. Numerical results presented in [13], 
indicate that for a given level of fidelity, or approxima- 
tion, the required times grow linearly with N, the num- 
ner of sites in the chain minus one. (We shall comment 
on these observations in Section 5, using exact results 
on almost perfect return [14] as applied to the uniform 
XX chain.) Thus Godsil [12] expressed the view that 
pretty good state transfer or APST will not be a satis- 
factory substitute for perfect state transfer in practice. 
We here wish to stress in this connection that when con- 
sidering non-uniform couplings, there are situations, as 
shown by an aforementioned example in Section 5, where 
high fidelity can be achieved in finite time irrespectively 
of the chain length. This suggests that the class of spin 
chains with APST, a much larger one than the PST fam- 
ily, should not be so radically deemed of impractical use 
at this point. 



2. 1-EXCITATION STATE TRANSFER IN XX 
SPIN CHAINS 

We shall consider XX spin chains with nearest- 
neighbor interactions. Their Hamiltonians H are of the 
form 

N-l N 

H = g £ Jwi°l°f+i + <ri°l+i) + 2 £ + !)' 

i=o i=o 

(2.1) 



where J/ are the constants coupling the sites I — 1 and I 
and Bi are the strengths of the magnetic field at the sites 
I (I = 0,1, ... , N). The symbols of, of, erf stand for the 
Pauli matrices which act on the l-th spin. 
It is immediate to see that 

1 N 

[H,~ 2 £(of + l)]=0 
z=o 

which implies that the eigenstates of H split in subspaces 
labeled by the number of spins over the chain that are 
up. In order to characterize the chains with APST, it 
will suffice to restrict H to the subspace spanned by the 
states witch contain only one excitation. A natural basis 
for that subspace is given by the vectors 

\e n ) = (0,0,...,1,...,0), n = 0,l,2,...,N, 

where the only " 1" occupies the n-th position. The re- 
striction J of H to the 1-excitation subspace acts as fol- 
lows 

J\e n ) = J„+i|e n+ i) + B n \e n ) + J n |e„_i). (2.2) 
Note that 

Jo = Jn+i = (2.3) 

is assumed. 

Consider the polynomials Xn (x) obeying the recurrence 
relation 

Jn+lXn+l(x) + B n Xn(x) + JnXn-l(x) = XXn(x) (2.4) 

with 

X-i=0, xo = l- (2.5) 
They satisfy the orthogonality relations 

JV 

^2w s Xn(x s )Xm(x s ) = S nm , (2.6) 
s=0 

where w s are the discrete weights taken to satisfy 

JV 

$>. = !• (2.7) 

s=0 

In what follows we shall take the eigenvalues x s in in- 
creasing order 

x Q < x\ < x 2 < ■ ■ ■ < x N . (2.8) 

Let 

\x s ) = £ Vw~sXn{x s )\e n ) . (2.9) 

n=0 

It is easily seen that these vectors \s) are eigenstates of 
J with eigenvalues x s 

J\x a ) =x s \x s ). (2.10) 
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Since both bases {\s)} and {|e„)} are orthonormal 
we also have 

N 



\ e n) =^2V™~sXn{x.s)\x S )- 



(2.11) 



s=0 



Let Pjv+i(a;) be the characteristic polynomial of J 
P N+ i(x) = (x- x )(x - xi) . . . (x - x N ). (2.12) 
The discrete weights are expressed as [10] 
h N 



w s = 



Pn{x s )P' n+1 {x s )'' 



s = 0,l,...,N (2.13) 



with ft, 7v = J\ J% ■ ■ ■ J% an( i Pn{x) = h^XNix). 



3. NECESSARY AND SUFFICIENT 
CONDITIONS FOR ALMOST PERFECT 
QUANTUM STATE TRANSFER 

By "almost perfect state transfer" (APST) we mean 
the following. 

Assume that the spin chain is prepared at time t = 
in the pure state |eo). This means that that at site n = 
the spin is up while all other spins are down. 

For arbitrary times t, this state will evolve into the 
state 



|e (t)> = e- lU \e ) 



(3.14) 



We demand that for any small e > 0, there exists a value 
of time t such that 



M*)> 



J4>(t) 



ejv) 



< c, 



(3.15) 



where (f>(t) is a real parameter which can depend on t. 
This means that the state \e (t)) can be as close to the 
state |ejv) as desired and that \eo) has thus undergone at 
time t an almost perfect transfer. 

The notation \t] — £ | 2 for two vectors £ and 77 stands as 
usual for 



l£-r?| 2 = £|a-%| 



(3.16) 



fe=0 



where are the expansion coefficients of the vector £ 
over a basis, say |efe): 



N 



(3.17) 



k=0 



Recall that the condition for 1-excitation PST reads 



M*)> 



-%4>(t) 



(3.18) 



this means that there then exists a time t for which the 
state |eo(i)) coincides (up to a phase factor e 1 ^) with 
the state |ejv)- 

In the case of almost perfect state transfer, there is no 
time t for which condition (3.18) is verified. Neverthe- 
less, it is possible to approach the state |ejv) with any 
prescribed degree of accuracy. From a practical point of 
view, there is no essential difference between perfect and 
almost perfect state transfer, owing to inevitable tech- 
nological and measurement errors. However, the APST 
conditions are much weaker than the PST ones. APST 
therefore widens the possibilities for constructing efficient 
quantum wires. 

Let us first derive the necessary condition for APST. 

Taking into account expansion (2.11), we have 



N 



N 



\x S ), 



s=0 

and hence 



l e W> = X y/uhXN(x s )\x s ) 

(3.19) 



s=0 



e 

N 

£* 

s=0 



* (t) |eo(t)>-M 



Xn(x s ) 



(3.20) 



It is easily seen that in order to fulfill condition (3.15) we 
need 



\xn(x s )\ = 1. 



(3.21) 



To convince oneself of this fact, suppose that (3.21) does 
not hold and assume that for some s = 0,1, ... ,N, we 
have xn{x s ) — a with \a\ 7^ 1. In this case, using the 
reverse triangle inequality we have 



<*<*>-«*• - X n(x s ) >(|a|-l) 2 >0 



and as a consequence the r.h.s. of (3.15) cannot be made 
arbitrarily small. 

Since xn(x) has only real coefficients (3.21) implies 
that xn{x s ) = ±1- From general properties of orthog- 
onal polynomials [10], [23] it follows, using (2.12), that 
(3.21) amounts to 



Xn{x s ) = (-1) 



N+s 



(3.22) 



As shown in [23], (3.22) implies that the matrix J is 
persymmetric, or mirror-symmetric i.e. 



RJ R — J, 

where R is the reflection matrix 



(3.23) 



R = 











Vi 



/ 
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We have thus obtained a necessary condition for APST 
namely, that the Jacobi matrix J corresponding to the 
spin XX Hamiltonian should be mirror symmetric. This 
coincides with one of the necessary conditions for PST 
[17]. The other condition for PST requires that 

x s+1 -x s = jM s , (3.24) 

where M s are arbitrary positive odd numbers. 

We shall now obtain the conditions on the spectrum of 
J for APST. 

The matrix J is Hcrmitian and mirror symmetric with 
Ji > and hence all its eigenvalues real and dis- 

tinct. 

Let fon(t) be the amplitude for finding the system in 
the state \n) at time f if was in the state |eo) at time 
f = 0, 

kn{t) = (e n \e- iJt \e ). (3.25) 
It is easily seen that [23] 

N 

f 0N (t) = ^- iXat XN(x s ). (3.26) 
Taking into account (3.22) we have 

N 

foN(t) = J2 w * e ~ iXat (- 1 ) N+S - ( 3 - 27 ) 

Mindful of the normalization (2.7), it is seen that con- 
dition | /oat (t) | « 1 is equivalent to the condition that 

e -ix a t^_^N+s _ & i<t> ( 3 28 ) 

for a fixed value t, where ss means "approximately equal 
with any prescribed accuracy". To paraphrase (3.28), 
there thus should be a value f of time for which the l.h.s. 
is as close as desired to a phase independent of s. 

This means that |/ojv(£)I is an almost-periodic func- 
tion. 

Recall that any almost periodic function /(f) is a for- 
mal trigonometric series [2] 

oo 

/(f) = ]T one*""*, (3.29) 

n— — oo 

where u> n are real parameters. For periodic functions 
f(t + T) = /(f), one has uj n = an d (3.29) is the 
ordinary Fourier series. For almost periodic functions 
there exist the so-called almost periods. This means that 
for every e > there exists a real parameter T = T(e) 
such that the inequality 

|/(f + T)-/(f)| <e (3.30) 

holds for all f. 

In turn, condition (3.28) is tantamount to an inequality 
involving the exponents that can be stated as follows: 



For every S > there exist a real parameter (j> and a 
value of time f such that 

\x s t - its + <j>\ < 6 (mod27r). (3.31) 

In more details, (3.31) can be rewritten in the form 

-6 < x s t - 7ts + <j> + 2ttM s < 8, (3.32) 

where M s are integers which may depend on s. 

This obviously amounts to a condition on the spectrum 
of J for APST to occur. The question is: what are the 
properties that the eigenvalues x s must possess to ensure 
that it is possible to find a time f and integers M s so 
that (3.32) is satisfied and hence APST realized. In other 
words, given the set of real numbers a s — cf> — its, this is 
asking: when is it possible to find values of f for which 
x s t is approximated in terms of integers by a s — 2irM s 
with any prescribed accuracy? 

The solution of this Diophantine approximation prob- 
lem (3.32) is given remarkably by the Kronccker theorem 
[15], [19]. 

In order to state it, let us introduce the following def- 
inition [15], [19]: A set on, i — 1, 2, . . . of real numbers is 
called linearly independent if for any n the only rational 
values of r\ , . . . , r n satisfying 

nai + r 2 a 2 H h r n a n = 

are r\ = r 2 = ■ ■ ■ = r n = 0. 

The Kronccker theorem can be formulated in two ver- 
sions, one a special case of the other. The first version is 
attributed to Kronecker himself. 

Version 1 Assume that the real numbers x s , s = 
0,1,... N are linearly independent over the field of ra- 
tional numbers. Let a , a\, . . . , ajy be fixed arbitrary real 
numbers. The Kronecker theorem states that for every 
S > there exist a real t and integers Mo, M\, . . . Mn 
such that the inequalities 

\x s t - a s - 2wM s \ <5, s = 0,1,2,..., N (3.33) 

hold. 

We see how this theorem directly applies to the APST 
problem. Indeed, provided that the eigenvalues are lin- 
early independent over the rational numbers, we may 
conclude from the Kronecker theorem that for every pa- 
rameter (f> it is possible to find a time f so that |/ojv(t)| 
is as close to 1 as desired. 

In many cases however, the eigenvalues not lin- 

early independent. This means that there can exist L 
independent relations of the type 

r^x + rfxi + ■■■ + r ( $x N = 0, i = 1, 2, . . . , L < N 

(3.34) 

(i) 

where arc integers such that for every i = 1, 2, . . . , L 
at least one of them is nonzero. (The use of integers is 
equivalent to that of rationals and more practical.) 

When additional relations such as (3.34) are present, 
the Kronecker theorem can be formulated as follows [19]. 
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Version 2 Assume that the real parameters x s , s — 
0,1,..., N are all distinct and moreover that there 
are L relations of the type (3.34) with nontrivial sets 

{r^ , . . . rjy } of integers. 

Then, the approximation condition (3.33) holds for ev- 
ery S > if and only if the real quantities at satisfy the 
conditions 

ao+r^ a±-\ \~ r $ a N = (mod2n), i = l,2,...,L 

(3.35) 

with the same integers rj 1 ' as in (3.34). 

For the proof of this statement see [19]. 

Using this (generalized) version of the Kronecker the- 
orem we can formulate the necessary and sufficient con- 
ditions for APST. 

General result. Let xq, x\, . . . , xjy be N + 1 distinct 
eigenvalues of the Jacobi matrix J corresponding to the 
XX spin chain (2.1). Assume that there are L < N rela- 
tions of the type (3.34) with nonzero integer parameters 

-(0 

' s • 

Then, the following conditions are necessary and suf- 
ficient for APST: 

(i) the Jacobi matrix J is mirror- symmetric, i.e. 

B s = B N _ S , J s = J N+1 _ S , s = 0,l,...N; (3.36) 

(ii) the L linear relations 

N 

^rf^Trs -</>) = (mod2Tr), i — 1,2, ... L (3.37) 

s=0 

are compatible. 

An obvious special case of this statement occurs when 
all eigenvalues x s are linearly independent over the field 
of rational numbers. This means that L = 0, i.e. that 
there are no additional relations such as (3.37) and that 
version 1 of the Kronecker theorem suffices to conclude to 
the occurrence of APST . In such instances , we see that 
the mirror symmetry of J and the linear independence 
of the eigenvalues represent the necessary and sufficient 
conditions for APST. 

Let us now compare the conditions for APST with the 
conditions for PST. It is known that the PST conditions 
are [17]: 

(i) the Jacobi matrix J is mirror-symmetric; 

(ii) Assuming that the eigenvalues x s are ordered so 
that xo < x\ < ■ ■ ■ < xn, the differences A s = x s+1 — x s 
are proportional to odd numbers: 

A s = K (2j s + l), (3.38) 

where j s are integers and re is a constant not depending 
on s. 

We see that mirror-symmetry is necessary for both 
PST and APST. This implies that the matrix J and 
hence the Hamiltonian can be uniquely reconstructed 
from the eigenvalues xq, x\, . . . Xn when either process 



occurs. Indeed, since (3.22) (which implies mirror- 
symmetry) specifies xn{x) at N + 1 distinct points, this 
polynomial can be obtained by Lagrange interpolation. 
Two orthogonal polynomials associated to the 3-diagonal 
matrix J are therefore known: xn(x) and the charac- 
teristic polynomial IlfLo( a ' — Xs ) °^ degree N + 1. As 
shown in [23], all the orthogonal polynomials can then 
be found by the Euclidean algorithm thereby providing 
their 3-term recurrence relation explicitly, i.e. giving all 
the couplings Ji and magnetic fields E>i. 

The difference in the conditions for PST and APST 
arises only in the restrictions on the spectrum or the 
eigenvalues of J. As is observed, these conditions are 
much more stringent for PST than for APST. 



4. SPECTRAL SURGERY AND APST 

A spectral surgery procedure was proposed in [23], as 
a way to construct XX spin chains with PST from such 
systems already known to possess the PST property 

Given the initial spectrum Xn = {xq, x\, . . . ,Xn}, un- 
der spectral surgery, one or several levels x^ , Xi 2 , . . . , Xi j 
are removed from the set A/y. 

The reduced set is Ajv-j = Xn \ Rj, where Rj = 
{x i± , Xi 2 , . . . , x^ } is the set of levels that have been 
removed. The spectral levels determine the (mirror- 
symmetric) Jacobi matrix J (and hence the Hamiltonian 
of the XX spin chain) uniquely. Thus the new spectral 
set Xft-j generates a new Jacobi matrix J of dimen- 
sion (N + 1 — j) x (N + 1 — j). Certain restrictions 
need to be imposed on the possible choices of sets Rj. 
Namely, Rj can contain any number of levels starting 

from the first one, say R^\ = { x o, x ii ■ ■ ■ i #Mi-i} or 
similarly, any number of levels from the last level, say 
R^Mz = { x n,xn-i, . . . ,x n -M 2 +i}- The only restriction 
is that there should be no gaps in the above sequences. 
Apart from these two elementary surgeries, levels can 
also be extracted from the middle of the spectral set 
Xn; in this case the requirement is the following: all 
such level sets should consist in the union of subsets 
R^ = {xi, Xi+i, . . . , iEi+L-i} of even length L and with- 
out gaps within R^ . Equivalently, one can say that only 
pairs Xi, Xi+i of neighbor levels can be removed from the 
middle of the set X^. 

Remarkably, the Jacobi matrix J corresponding to a 
surgered set X^-j, is obtained from the initial Jacobi 
matrix by j Christoffel transforms (see [23] for details). 

That is, after j Christoffel transforms, we obtain a new 
Jacobi matrix J which is mirror-symmetric and satisfy 
the PST property. This observation allows to generate 
new explicit examples with PST without the need to per- 
form the inverse spectral problem algorithm. This is ad- 
vantageous since the formulas of the Christoffel transform 
are rather simple. 

One such example was already given in [23]. Let 
N be odd. We start with the uniform grid Xn = 
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{-N/2, -N/2 + 1, . . . , N/2} and then remove 2L levels 
symmetrically from the middle of the set X^ 

X N - 2L = {-N/2, -N/2 + 1, L — 3/2, —L - 1/2, 

L + 1/2, L + 3/2, . . . , N/2 - 1, N/2}. (4.1) 

Such a spectral set corresponds to the model of the XX 
spin chain with PST proposed in [21]. 

Consider what is the effect of spectral surgery with 
respect to APST. 

It is almost obvious that under an arbitrary admissible 
surgery procedure the APST property survives. Indeed, 
if all eigenvalues Xi were linearly independent over the 
field of rational numbers then any reduced set X^-j will 
obviously satisfy the same property. 

Assume that there exists a set of linear relations (3.34) 
on the spectral levels x s . APST is possible iff the set of 
linear relation (3.37) is compatible (it is assumed a priori 
that the matrix J is mirror-symmetric). 

Under spectral surgery some levels of the set Ajy are 
removed. This means that some of relations (3.34) will 
disappear (or remain the same). This implies that the 
new Jacobi matrix J will also possess APST. 

Interestingly, this method can sometimes generate a 
chain with APST, even if the initial chain lacks this prop- 
erty. Indeed, the absence of APST for a given mirror- 
symmetric XX chain means that relations (3.37) are in- 
compatible. When we remove some of the levels the 
number of relations (3.34) can be reduced which in prin- 
ciple can lead to the compatibility of the reduced set of 
relations (3.37). Of course one should seek to realize this 
in concrete examples. We plan to investigate this prob- 
lem in the future. 



5. EXAMPLES AND SPECIAL CASES 

We shall present a number of examples and special 
cases in this section: 

i) we shall indicate how the central result specializes 
in the absence of magnetic fields; 

ii) we shall record the circumstances for APST in the 
uniform XX chain and shall discuss the relation between 
the waiting time and the length of the chain for this par- 
ticular system, using results on almost perfect return (to 
the point of origin); 

iii) we shall show, not too surprisingly, that the con- 
ditions for PST are a special case of the conditions for 
APST; 

iv) we shall consider in details the case N — 4 corre- 
sponding to spin chains with 5 sites. In this instance, 
it is possible to give a very explicit analysis of APST 
and to estimate waiting times using standard methods of 
Diophantine approximation. The conditions for all the 
possible cases (PST, APST and neither) are derived. 

v) we shall introduce a remarkable model based on 
the para-Krawtchouk polynomials that can exhibit PST, 
APST or neither depending on the value of one of its 



parameters ; it will furthermore be seen that the wait- 
ing times can be estimated and found to be finite inde- 
pendently of the size of the chain for arbitrarily good 
fidelities in some special APST situations. 

i. Absence of magnetic fields. 

Consider the special situation corresponding to zero 
magnetic fields B s = 0. In this case the Jacobi matrix J 
has only two nonzero diagonals. The corresponding or- 
thogonal polynomials are symmetric P n (—x) = —P n (x) 
and hence the eigenvalues satisfy the properties 



x s + x N - s = 0, s = 0, 1, 



,A 



(5.1) 



where it is assumed that xq < x\ < . . . xn- 

Consider first the case N odd. Then one has (N + 
l)/2 independent relations (5.1) for s = 0, 1, . . . , (N — 
l)/2. They coincide with relations (3.34) where only two 
integers are nonzero and given by r s = rjv- s = 1- 
Relations (3.37) then reduce to the single condition 



■wN - 2(f> = (mod 2tt) 



(5.2) 



from which we find 



-tt/2 if A = 4m + 3 
tt/2 if A = 4m + 1 



(mod 2tt). (5.3) 



If one assumes that the eigenvalues Xo,x\, . . . , £(jv-i)/2 
are linearly independent, then the Kronecker theorem 
guarantees with (f> — ±7r/2 that the corresponding spin 
chain will exhibit APST. 

If N is even, we shall have the N/2 independent rela- 
tions (5.1) for i = 0, 1, . . . , N/2 — 1 and in addition 



x N/2 



0. 



Conditions (3.37) are compatible iff 







if N/2 
if N/2 



is even 
is odd 



(mod 271"). 



(5.4) 



(5.5) 



Again, if the eigenvalues Xo, x\, . . . x N / 2 -i are linearly 
independent, the APST property is present. 

ii. The uniform XX chain with zero magnetic 
fields. 

The APST (or pretty good state transfer) of the XX 
chain with uniform couplings and no magnetic fields, i.e. 
J s = 1 and Bs — for all s, have recently been sorted 
out in [13] . It has been found there that a chain of length 
N+l admits APST if and only if N = p-2 or 2p-2 where 
p is prime or iff N = 2 m — 2. As a matter of fact, the proof 
of that result given in [13], effectively proceeds from the 
application of version 2 of the Kronecker theorem and 
the use trigonometric relations between the eigenvalues 

x s = 2 cos(tts/(A + 2)), s = 1, 2, . . . , N + 1 (5.6) 

which are given by the roots of the Chebyshev polyno- 
mials of second kind for this model. 

In the same reference [13] (see also [12]), it is pointed 
out, on the basis of numerical analysis, that the waiting 
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times for APST will grow with the size of the chain. This 
has cast doubts on the practical use of APST. Our next 
examples will hopefully dissipate this negative view by 
showing that there are models exhibiting APST where 
the waiting times stay finite as the size is grown. 

Before we leave this uniform XX chain, we would like 
to offer the following comment on the waiting times as 
the chain is taken to have very large length. 

The quantity 

/oo(t) - (eo|e^ tJ |eo) (5.7) 

represents the return amplitude, i.e. the amplitude to 
return to the initial (input) state after a time t. Since 

|e-* tJ |e > - e**jejv>| = \e itJ \e N ) - e"^|e )| , 

if we have APST from |eo) to |ejv), we also have APST 
from |ejv) to |eo). As a consequence, in a chain with 
APST, we must observe that there is almost perfect re- 
turn, that is that there are times t n , n = 0, 1, 2, . . . such 
that 

lim |/oo(tn)| = l. (5.8) 

n— >oo 

The conditions for almost perfect return have been an- 
alyzed in [14] and can be simply stated: An XX spin 
chain will exhibit almost perfect return if and only if 
its one-excitation Hamiltonian J has a pure point spec- 
trum. As is readily seen, the spectrum of the uniform 
XX chain becomes continuous as N , the number of sites 
minus one, becomes infinite (see [14] for more details.) It 
follows that almost perfect return does not occur in the 
case of a semi-infinite Heisenberg XX chain and hence 
not surprisingly an infinite time is required for APST in 
this model. 

iii. PST as a special case of APST. 

If a set {xq, Xi, . . . xn} of eigenvalues satisfy the APST 
conditions, it is clear that the affine transformed eigen- 
values x s = ax s + (3 will also satisfy the APST condi- 
tions. In particular, we can always assume that xq = 
and x\ = 1. The PST condition (3.38) then implies that 
x 2s = K 2s are even integers while X2 S +i = K 2s +i are odd 
integers. We have N linear relations 

K s x x =x s , s = 0,2,3,..., N. (5.9) 

With 4> = we see that the conditions (3.37) become 

K s = s(mod 2). (5.10) 

This is equivalent to the condition that the parity of the 
integers K s coincides with the parity of s, i.e. with the 
PST condition. We thus see that the PST condition is a 
special case of the APST condition, as expected. 

A note is in order here. In general the translation 
x s — > x s + (3 of the energy spectrum has no physical 
meaning: it is always possible to redefine the value of 
the ground state energy and for instance to put it equal 
to zero as noted. (Similarly the scale can be arbitrarily 



chosen using the dilation factor a.) This freedom relates 
to the values that the phase (f> takes modulo 2ir. For a 
given spectrum, the APST conditions might restrict (f>, 
however shifts in the ground state energy will allow <j> to 
take arbitrary values. 

iv. Spin chains with 5 sites and without mag- 
netic field 

Consider the the XX spin chain without magnetic 
fields (Bo = B\ = ■ ■ ■ = Bjy = 0) that contains 5 sites. 
In this case it is possible to express the APST conditions 
explicitly and to describe 3 possible scenarii: PST, APST 
and neither. 

This corresponds to N = 4. There are only two dis- 
tinct exchange constants: J\ and J 2 . Indeed, due to the 
mirror-symmetry condition (3.36) we have (for N = 4) 
that Ji — J\ 1 J3 = J 2 . 

The spectrum of the corresponding Jacobi matrix J 
can easily be found: 

x = —a, X\ = —b, x 2 = 0, x 3 — b, x 4 = a, (5.11) 

where a = \/ J 2 + 2 Jf , b = J\ are two positive param- 
eters. Clearly a > b so the eigenvalues x s in (5.11) are 
given in increasing order. 

This is a special case of the situation considered in i. 
The sufficient condition for APST is that the parameters 
a and b be linearly independent over the rationals. This 
is equivalent to the statement that the ratio a/6 is an 
irrational number. We thus see that if \J J\ + 2 J|/ J\ is 
an irrational number, then the spin chain with 5 nodes 
has the APST property. 

Suppose now that the ratio a/b = n/m is a rational 
number where n, m are co-prime integers. This means 
that there are two conditions on the eigenvalues 

mxo — nx\ = mx4 — nx 3 = (5.12) 

Then the relations (3.37) which are required for APST 
(with <f> = 0) reduce to the only condition that n be even 
and m be odd. 

But this is equivalent to the PST condition. Indeed, 
we have 

^ = - - 1 = M/m, (5.13) 
b m 

where M and m are both odd. Observe that the dif- 
ferences between the eigenvalues are Ai = x\ — xq = 
a - b, A 2 = x 2 - x\ = b, A 3 = x 3 - x 2 = b, A 4 = 
X4 — X3 = a — b. Hence (5.13) coincides with the PST 
condition for a chain with 5 sites. 

In all other cases, i.e. when n is odd and m is odd, or 
when both n and m are odd, the APST condition is not 
valid and the XX spin chain does demonstrate neither 
PST nor APST. 

In the APST situation (i.e. when a/b is irrational), it 
is interesting to estimate the waiting times. It is easy to 
calculate the amplitude /ojv(t) directly from (3.27): 

a 2 -b 2 b 2 1 
/oiv(i) = — !~ 2^2 cosa * ~ 2 C ° S& * ^ 5 ' 14 ' ) 
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This is a simple example of almost periodic functions: 
two harmonic oscillations with frequencies uj\ — a and 
uj 2 = b do not constitute a pure periodic function. But 
the Kronecker theorem guarantees that there exists a se- 
quence t n such that foN(t n ) — >■ 1 when n —¥ oo. 

There are several simple algorithms to estimate the 
quantities t n . 

One of them proceeds through the expansion of a/b 
into a continued fraction [15] 



a/b = Co + 



1 



(5.15) 



Ci 



C2 + 



where the quotients Co>Ci>--- are integers. The corre- 
sponding convergents p n /q n of the continued fraction are 
rational numbers which provide the best approximation 
of the irrational number a/b [15]. In our case this method 
works well if all (or at least infinitely many) convergents 
have the property that p n is even and q n is odd. Oth- 
erwise one can apply the method of Farey fractions [15] 
or, equivalently, the method of intermediate convergents 

[18]- 

This method can described as follows. Starting with 

the convergents {po/Qo,Pi/Qi, ■ ■ ■ } wc compute the so- 
called mediants (or intermediate convergents) by the for- 
mula [18] 



pn _ Pn+1 ±Pn 
<7n q n +i ± Qn 



(5.16) 



The convergents p n /q n together with the intermediate 
convergents p n /q n form a set of generalized convergents 
w n = u n /v n , n = 0, 1, 2, . . . , where u n = p n or u n = p~ n 
(similarly v n = q n or »„ = q n ). From the elementary 
properties of the convergents it follows that all fractions 
u n /v n are simple (i.e. u n and v n are co-prime). More- 
over, it is possible to show [18] that the generalized con- 
vergents provide the general solution to the best approx- 
imation problem: they satisfy the inequality 



\z - u n /v n \ < v~ 



(5.17) 



It is easy to show that there are infinitely many general- 
ized convergents u n /v n with the desired property: u n is 
even while v n is odd. 

Indeed, assume that there is only a finite number of 
the convergents p n /q n with the desired property. This 
means that starting with some n = M the numerators 
p n are all odd; the corresponding denominators q n may 
be either even or odd. Necessarily, for any neighbor pair 
Pn/qn and p n+1 /q n+ i, the denominators q n and q n+ i 
should be of opposite parity. Indeed, if one assumes 
that q n and q n +i have the same parity then their me- 
diant (p n +p n +i)/(q n + q-n+i) will be reducible (with 
both numerator and denominator even) which is impos- 
sible. Thus the mediants p n /q n will have the desired 
property, because p n = p n +i ± p n is even as a sum of 
two even numbers and q n = q n +i ± q n is odd as a sum of 



two number with opposite parities. This means that for 
any irrational number z = a/b there exists an infinite se- 
quence of generalized convergents c„ = u n /v n such that 
linin^oo c„ — > z and u n is even while v n is odd. More- 
over, the best approximation property (5.17) holds for all 
such convergents. 

Given this sequence, let us put 



TTVn 

b 



(5.18) 



The amplitude /ojv(in) then becomes 

<r Ir Ir / , M I 



, . a" — o o ( a\ 1 

f0N(tn) - + ^ COS + -. (5.19) 



This can easily be simplified to 

b 2 



/<w(t n ) = 1 - — sin 2 (^„e„/2), (5.20) 



where 



- n 



a Un 
b V n 



(5.21) 



is the accuracy of the rational approximation of the irra- 
tional number a/b. 

By property (5.17) we have that 



foN ~ 1 



ir 2 b 2 



(5.22) 



It is seen from this formula that /oiv(in) converges to 1 
when n — > oo. Moreover, this formula allows to estimate 
the accuracy of this approximation if the (generalized) 
convergents u n /v n are known explicitly. 

It is seen that /oat (t n ) — > 1 when t n — > oo . This means 
that the sequence {t n } is associated to APST. 

Consider the special case J\ = J2 = 1 so that all cou- 
plings are now equal. For the homogeneous Heisenberg 
chain (i.e. with Ji = 1 for all i = 0,1,2, N) it is 
known [11] that PST only occurs when there are 2 or 3 
sites. We are thus in a situation of APST for a uniform 
chain with 5 nodes. 

In this case a = 6 = 1. The continued fraction 
representation is 



V3= 1 + 



1 



(5.23) 



1 + 



2 + 



1 + ... 



The first convergents are 



Pn = f , 5 7 19 26 71 
q n \ ' ' 3' 4' 11' 15' 41'' 



(5.24) 



It is seen that the first two convergents with the desired 
property (i.e. an even numerator and an odd denomi- 
nator) are 2/1 and 26/15. Computing the mediants ac- 
cording to (5.16), we obtain the additional intermediate 
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convergent 12/7. Hence, the first 3 desired (intermedi- 
ate) convergents are 



Un 
V„ 



2 12 26 

I'T' 15 



(5.25) 



The first fraction, 2/1, yields an approximation to the 
difference between |/ojv(i)| and 1 with an accuracy of 
3- 10 _1 ; the second fraction 12/7 yields an approximation 
with accuracy 2 • 10~ 2 . 

So, already in the case of a spin chain with 5 sites we 
see that all 3 possibilities: PST, APST or neither can 
occur. APST is generic, i.e. APST happens for almost 
all possible values of J\ and J 2 . 

We shall consider next a XX spin chain with an ar- 
bitrary (even) number of spins which demonstrates the 
similar property. 

v. Spin chain corresponding to the para- 
Krawtchouk polynomials. 

Assume that N is odd and that the spectrum x s is 
such that (up to an affine transformation) 

x 2s = 2s, x 2s+1 = 2.S+7, s = 0, 1, . . . , (iV-l)/2, (5.26) 

where 7 is a real parameter such that < 7 < 2. 

This situation can be realized in spin chains that are 
obtained from the para-Krawtchouk polynomials with 
the recurrence coefficients [24] 



B n — 



7 + AT- 1 



2 _ n(N + 1 - n)((N + 1 - 2n) 2 - 7 2 ) 



4(N -2n)(N ~2n + 2) 



(5.27) 



The corresponding Jacobi matrix is mirror symmetric. 
When 7=1 the grid is uniform and the Jacobi matrix J 
which is associated to the ordinary Krawtchouk polyno- 
mials generates PST [1]. 

It is clear that all eigenvalues can be expressed as linear 
combinations of x\ and x 2 . Indeed, 

x 2s = sx 2 , x 2s+1 = x x +sx 2 , s = 0, 1,2, . . . (N- l)/2 

(5.28) 

Hence, according to our central result (3.37), in order to 
have APST, the following conditions must be realized: 



a 2s 



sa 2 , 



a 2s+i = a i + sa 2, (mod 2ir) (5.29) 



where a s = 7rs — <fi. These relations are compatible iff 
= (mod2vr). 

If the parameter 7 is irrational then there are no addi- 
tional relations for the eigenvalues x\, x 2 and thus pro- 
vided (5.29) is verified, APST will occur. 

Equivalently, this means that there exists a sequence 
{t n , n — 1,2,...} such that 



lim |/o;v(in)| = 1 



(5.30) 



As in the previous case it is possible to determine the 
times t n and estimate the corresponding rate of conver- 
gence in (5.30) using standard Diophantine approxima- 
tion methods. . 



Indeed, we can always choose an infinite set of (inter- 
mediate) convergents {ui/i>i, u 2 /v 2 , . . . , u n /v n , . . . } with 
the property that the numerators u n are all odd. Intro- 
duce the value 



7 



(5.31) 



which describes the accuracy of the Diophantine approx- 
imation of the parameter 7. We already saw that 



Choose 



t n = TTV n . 



(5.32) 



(5.33) 



Substituting this in (3.27) and taking into account that 
u n is odd we have 



(JV-l)/2 (JV-l)/2 
f0N(t n ) = ~ ]T W 2s + 

(JV-l)/2 (JV-l)/2 



W 2s+ 1 



J2 W2s 



J2 W2 > 



(5.34) 



s=0 



s=0 



In [24] it was shown that for the para-Krawtchouk poly- 
nomials the identities 



(JV-l)/2 (JV-l)/2 

^2 w 2s = ^2 w 2s+1 = 1/2 



s=0 

hold. Hence 



s=0 



foN(tn) 



(1 + e -<™»e») . 



Now 



\f0N(tn)\ = C0s(7TO„£„/2) W 1 - — 



(5.35) 



(5.36) 



(5.37) 



where we have used (5.32). 

Formula (5.37) gives a good approximation of the am- 
plitude |/o at I if v n is sufficiently large (depending on the 
physical requirements on the accuracy). E.g. for v n > 10 
we get an accuracy of 1%. 

Consider, e.g. the value 7 = 

From (5.23) we find the first appropriate (i.e. with all 
numerators u n odd) convergents 



u n _ J 5 7 19 45 71 
~ 1 ' 3' 4' IT '26' 41 : 



-}■ 



(5.38) 



Already the third convergent 7/4 gives an accuracy of 
about 1%. This is because the actual accuracy £3 w 
0.0179 of the third convergent is smaller than what the 
rhs of formula (5.32) gives. Hence the first waiting time 
corresponding to an accuracy of 1% in the amplitude is 
h = 4?r. 
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Consider now the case where 7 ~ p/q is a rational 
number, with p and q coprime integers. In this case there 
is an additional linear relation between x\ and x 2 : 

2qx\ — px2 (5.39) 

and hence a similar relation should hold for a\ and a 2 - 

2qai = pa 2 (mod 2n) (5.40) 

where a\ = it, a 2 — 2tt. If the numerator p is odd then 
(5.40) holds for any pair of coprime integers p,q and 
hence the APST condition is satisfied. In fact we are 
then in the situation where not only APST but even PST 
occurs. Indeed, this corresponds to the model with PST 
that we derived and discussed in [24]. 

If the numerator p = 2j is even for some integer j 
(in this case, necessarily the denominator q is odd) then 
relation (5.39) becomes qx\ — jx 2 and hence (5.40) reads 

q-rr - 2jir = (mod 2tt). (5.41) 

Obviously, (5.41) cannot hold when q is odd and there is 
no PST nor APST in this case. 

In summary the picture is as follows: 

(i) when 7 is an irrational number, APST is observed; 

(ii) when 7 = p/q is rational with p odd there is PST; 
(hi) when 7 = p/q is rational with p even neither PST 

nor APST happens. 

Let us make the following remark in connection with 
the last two examples. As mentioned in subsection iii, in 
the discussion of the uniform XX chain, almost perfect 
return occurs in chains of the type we have been con- 
sidering, this is so whenever the spectrum of the Jacobi 
matrix J is discrete. Clearly, this will be the case for 
any XX spin chain with nearest-neighbor interactions 
that has a finite number of sites. It should be stressed 
that finite length does not similarly imply APST; as we 
have seen, there are examples of systems with finite 1- 
excitation Hamiltonian J that do not possess APST. 

6. CONCLUSIONS 

Let us recapitulate the essential elements of our anal- 
ysis of almost perfect state transfer in spin chains. It 
builds on the knowledge that XX spin chains with prop- 
erly engineered couplings and magnetic fields can effect 
the transport of states from one end to the other with 
probability 1 over certain times - these are chains with 
perfect state transfer (PST). In view of the fact that there 
are always manufacturing or measurement imprecisions, 
our study aimed to characterize the models where al- 
though not perfect, state transfer could be realized with 
a probability very close to 1 - such chains have been said 
to show almost perfect state transfer (APST). 

We have thus undertaken to categorize all XX spin 
chains with arbitrary nearest-neighbor couplings and 
magnetic fields that would exhibit APST. The Kronecker 



theorem in Diophantine approximation proved essential 
in this investigation. The necessary and sufficient condi- 
tions that have been found for a chain to admit APST 
bear on the 1-excitation restriction J (a tri-diagonal ma- 
trix) of the chain Hamiltonian. Not unlike what is needed 
for PST, the requirements for APST entail a symme- 
try condition on J as well as spectral restrictions. As it 
turns out, like for PST, J must be mirror-symmetric that 
is, invariant against reflection with respect to its anti- 
diagonal. The conditions on the eigenvalues of J turn 
out (in general) to be much less demanding for APST 
than for PST. They amount roughly to the conditions 
for which the Kronecker theorem applies and imply that 
at least a subset of these eigenvalues be linearly indepen- 
dent over the field of the rational numbers. Spin chains 
with APST therefore much enlarge, in principle, the class 
of such systems that can be exploited as quantum wires 
since their fidelity can be as good as is technically rele- 
vant. 

It has also been demonstrated how chains with APST 
can be constructed from a chain that already admits (or 
possibly not) APST, by the removal of single excitation 
energy levels from the parent system via Christoffcl trans- 
forms. This offers a constructive method to broaden the 
catalog of systems with APST. 

A number of examples and special cases interesting in 
their own right have also been presented and discussed. 
They illustrate situations where PST or APST occur or 
where neither can happen. One of the models introduced 
have provided counter-examples to the view borne out in 
particular from the examination of the Heisenberg XX 
chain with homogeneous couplings, that APST requires 
excessively long times in extended wires. 

We trust this report provides interesting information 
on the transfer of state in XX spin chains with nearest- 
neighbor couplings and suggests that APST requires fur- 
ther analysis. 

Acknowledgments 

AZ thanks Centre de Recherches Mathematiques (Uni- 
versite de Montreal) for hospitality. The authors 
would like to thank M.Christandl, M.Derevyagin and 
A.Filippov for stimulating discussions. They also ac- 
knowledge communications from M.Bruderer, D. Bur- 
garth and S. Severini and are also very grateful to a 
referee for a very constructive and stimulating review. 
The research of LV is supported in part by a research 
grant from the Natural Sciences and Engineering Re- 
search Council (NSERC) of Canada. 



11 



[1] C.Albanese, M.Christandl, N.Datta, A.Ekert, Mirror in- 
version of quantum states in linear registers, Phys. Rev. 
Lett. 93 (2004), 230502. 

[2] A.S.Besicovitch Almost periodic functions, (Dover, 1954) 

[3] S. Bose, Quantum communication through spin chain dy- 
namics: an introductory overview, Contemp. Phys., 48, 
(2007), 13 - 30. 

[4] M. Bruderer, K. Franke, S. Ragg, W. Belzig, D. 
Obreschkow, Exploiting boundary states of imperfect spin 
chains for high-fidelity state transfer, Phys. Rev. A 85 
(2012), 022312. arXiv: 11 12.4503 

[5] D.Burgarth, Quantum State Transfer with Spin Chains, 
arXiv:0704.1309. 

[6] D.Burgarth S.Bose, Conclusive and arbitrarily perfect 
quantum state transfer using parallel spin chain chan- 
nels, Phys. Rev. A 71, (2005) 052315 . arXiv: quant um- 
ph/0406112v4. 

[7] C.Burrell and T. Osborne, Bounds on information prop- 
agation in disordered quantum spin chains, Phys. Rev. 
Lett. 99, (2007), 167201 ArXiv: quant-ph/0703209v3 
[8] C.Burrell, J.Eisert and T.Osborne, Information propa- 
gation through quantum chains with fluctuating disorder, 
Phys. Rev. A 80 (2009), 052319 . ArXiv: 0809.4833vl. 
[9] G. De Chiara, D. Rossini, S.Montenegro and R. Fazio, 
From perfect to fractal transmission in spin chains, Phys. 
Rev. A 72, 012323 (2005) arXiv: quant-ph/0502148v2. 

[10] T. Chihara, An Introduction to Orthogonal Polynomials, 
Gordon and Breach, NY, 1978. 

[11] M. Christandl, N.Datta,Tony C. Dorlas, A.Ekert, A.Kay 
and A. J. Landahl, Perfect transfer of arbitrary states in 
quantum spin networks, Phys. Rev. A 71 (2005), 032312. 
arXiv:quant-ph/0411020 

[12] C.Godsil, State Transfer on Graphs, Discrete Math. 
312(1): 129-147 (2012). arXiv: 1102.4898 



[13] C.Godsil, S.Kirkland, S.Severini, J.Smith, Number- 
theoretic nature of communication in quantum 
spin chains, Phys. Rev. Lett. 109 (2012), 050502; 
arXiv:1201.4822. 

[14] A.Griinbaum, L.Vinet and A.Zhedanov, Birth and death 
processes and quantum spin chains, arXiv: 1205.4689. 

[15] G.H.Hardy and E.M.Wright An introduction to the the- 
ory of numbers, 6th ed. Oxford University Press, 2008 

[16] A.Kay, Perfect State Transfer: Beyond Nearest-Neighbor 
Couplings, Phys. Rev. A 73 (2006), 032306 . ArXiv: 
quant-ph/0509065v2. 

[17] A.Kay, A Review of Perfect State Transfer and its Ap- 
plication as a Constructive Tool, Int. J. Quantum Inf. 8 
(2010), 641-676; arXiv:0903.4274. 

[18] . S.Lang, Introduction to Diophantine Approximations, 
Springer- Verlag, 1991. 

[19] B.M.Levitan and V.V.Zhikov, Almost periodic functions 
and differential equations, Cambridge University Press, 
1982. 

[20] C.Marletto, A.Kay and A.Ekert, How to counteract 
systematic errors in quantum state transfer, arXiv: 
1202.2978vl. 

[21] T. Shi, Y.Li , A. Song, C.P.Sun, Quantum-state trans- 
fer via the ferromagnetic chain in a spatially modulated 
field, Phys. Rev. A 71 (2005), 032309, 5 pages, quant- 
ph/0408152. 

[22] R.Ronke, T.Spiller, I.D'Amico, Long-range interactions 

and information transfer in spin chains, J. Phys.: Conf. 

Ser. 286 (2011), 012020. arXiv: 1101.4509vl 
[23] L.Vinet, A.Zhedanov, How to construct spin chains with 

perfect spin transfer, Phys.Rev. A 85 (2012), 012323. 
[24] L.Vinet and A.Zhedanov, Para-Krawtchouk polynomials 

on a bi-lattice and a quantum spin chain with perfect state 

transfer, arXiv:1110.6475v2. 



